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Abstract 

A simple model for the formation of a straight cosmic string, wiggly or un- 
perturbed is considered. The gravitational field of such string is computed in 
the linear approximation. The vacuum expectation value of the stress tensor 
of a massless scalar quantum field coupled to the string gravitational field is 
computed to the one loop order. Finally, the back-reaction effect on the gravi- 
tational field of the string is obtained by solving perturbatively the semiclassical 
Einstein's equations. 
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1 Introduction 



Cosmic strings are macroscopic topological defects that may have been produced at 
phase transitions in the early universe ^ 0]. They are predicted by gauge theories 
with spontaneous breakdown of symmetry whenever the unbroken subgroup contains 
a discrete symmetry as is the case of many grand unified theories (GUT), although 
not in the simplest unified group SU{5). The GUT strings would have been produced 
when the universe was 10~^^sec old and had a temperature of T ~ 10^^~^^GeV. The 
gravitational field of such strings may seed structure and, in fact, a network of these 
strings is an alternative to inflation for the generation of the universe structure ^, 

There are two types of gravitational quantum effects associated to cosmic strings 
(or to any other topological defect) which result from the interaction of the string's 
gravitational field with any quantum field {i.e. matter) present, namely, particle cre- 
ation and vacuum polarization. When a string forms, a sudden change in the gravi- 
tational field takes place which may translate into copious quantum pair production 
of particles in a way similar as electron-positron pairs are created by external elec- 
tromagnetic fields. This effect has been considered by several authors in different 
settings which go from various models of string formation ^ |^, ^ to oscillating 
string loops The main conclusion is that even though very energetic particles 



may be created the cosmological significance of these is very small compared to the 
background radiation at the time of formation. Typically, the ratio of the energy 
density of particles created by the formation of the strings and the energy density 
of the radiation is of order iV^(G/i)^, where N is the number of particle species, G 
Newton's gravitational constant, and n is the energy per unit length of the string; 
for GUT strings /i ~ lO^^g/cm, the square of the GUT mass, and thus Gfi ~ 10~^. 
For particles created by oscillating loops such ratio is much higher, N^^'^iGuY ^ but 
still cosmologically insignificant unless, of course, the number of particles is absurdly 
large; at GUT time we expect ~ 10^. 

Vacuum polarization effects due to quantum fluctuations of matter fields have 
been much less studied. This is due, in part, to the fact that one does not expect 
very significant changes in the classical string network picture as a consequence of 
such effects, but also in part because the computation of such effects is difficult 



T^ , p!3[ . The vacuum expectation value of the stress tensor of quantum fields around 
the string is generally different from zero even for a static string, consequently it is 
the source of a gravitational field which, in turn, modifies the classical gravitational 
field of the string. This is the back-reaction effect of quantum matter on the classical 
gravitational field. For example, the gravitational field outside an unperturbed {i.e. 
not wiggly) straight static string is described by the metric given by flat spacetime 
with a deficit angle in the plane perpendicular to the string W\ O, Il6|, O, |19 



The quantum stress tensor for conformally coupled scalar fields has been computed 
exactly by Linet [BO] and by Helliwell and Konkowski plj who found that the energy 
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density goes like Nh Gfi r~^, where r is the radial distance from the string axis. Such 
energy density creates in the weak field approximation a Newtonian potential outside 
the string of the order of $~ A^^ Gfir~'^. In fact, Hiscock [|1^ solved the semiclassical 



Einstein's equations to the linear order in this case and found the back reaction in the 
gravitational field of the cosmic string. The result is that the spacetime surrounding 
the string is no longer fiat with a deficit angle: the two- surface perpendicular to the 
string is an hyperboloid (rather than a cone) and the corrections to the fiat metric 
are of the order just described. Two consequences of this are clear, one is that a 
static string will exert Newtonian forces on surrounding non relativistic particles, the 
other is that when two cosmic strings approach they should feel increasingly strong 
attractive forces. The relative significance of these effects will be discussed later on. 

In this paper we compute the back reaction on the gravitational field of a cosmic 
string when it is formed. Now besides pure vacuum polarization effects, as in the 
case of a static string, we have an effect due to the creation of particles which also 
contributes to the quantum stress tensor. The computation is made within the linear 
approximation of the semiclassical correction to Einstein's equations, i. e. we assume 
that the spacetime metric departs from the Minkowski metric by linear terms h^^ 
only. This approximation is appropiate whenever the energy per unit length of the 
string satisfies that Gfi <^ 1, which is certainly true for GUT strings. Also the vacuum 
expectation value of the stress tensor for matter fields is computed perturbatively to 
first order in h^iy and to the one loop order using the results by Horowitz and 
Jordan P3[ when the background metric is fiat. More general results for the stress 



tensor in conformally fiat backgrounds are known |2^, |2^, ^ but here we ignore the 
effect due to the cosmological expansion. To this order the stress tensor does not 
include the energy of the particles created, which is second order in h^^, however 
particle creation effects such as those due to transition elements from vacuum to two 



particles states are included in this linear order, see the discussion in ref. . 

To implement string formation we use a very simple macroscopic model for the 
classical stress tensor of the string. The field dynamics of string formation is a too 
complicated issue to be described at microscopic level and one is restricted to consider 
rough macroscopic models for this process. The models are either implemented by 
giving the gravitational field of the string |^, ^, ||, ^ or by giving a prescribed form 



of the string stress tensor ||10| , |Tl|, |28|. The results for particle creation are rather 



independent of specific models and thus we take the simplest model based on the 



stress tensor given in ref. [0. In this model the string is always there, first as a dust 



straight line source and it is the string tension what grows in time, in this way the 
stress tensor is automatically conserved. 

The plan and a summary of the main results of the paper are the following. 
In section 2 we describe the model of string formation and derive its gravitational 
field in the linear approximation. We allow for the possibility that the string be a 
straight wiggly string. Straight wiggly strings are segments of long strings with small 
scale structure whose effective stress tensor may be described by a straight string 
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with effective mass per unit length larger than the unperturbed one {i.e. with no 
small scale structure) and an effective tension which is also less than the mass per 
unit length They appear as long strings in the numerical simulations of string 
evolution ||3y, ^ and may be the seeds of large scale structure [|, ^ As a 



consequence of linearity we can split the stress tensor in two parts, one which is 
static and whose gravitational field is easily solved and a time dependent part, whose 
gravitational field is found by solving the corresponding initial value problem. For 
simplicity we also use the sudden approximation, this introduces divergences in the 
gravitational field along the future light cone which must be regularized. To get rid of 
gauge effects which appear in the metric tensor we also compute the Riemann tensor. 
The evolution at large times of this tensor differs, of course, if the string is wiggly 
or unperturbed. In the first case we have Newtonian like potentials which help in 
building wakes |^ whereas in the second the curvature tends to zero corresponding 
to fiat space with a deficit angle. 

In section 3 the vacuum expectation value of the stress tensor outside the string 
for matter conformally coupled to the string's time dependent gravitational field is 
derived. It is seen that after a short transient period the stress tensor settles down 
to the values of a static, wiggly or unperturbed, string. In both cases, the energy 
density goes like NhGfir~'^ as one expects from dimensional arguments. Note that, 
for GUT strings, the string radius is tq ~ 10~'^°cm, corresponding to the Compton 
wavelength of the GUT mass and this radius set bounds on the energy density. 

In section 4 the semiclassical Einstein's equations for our problem are solved to 
find the back reaction of the quantum matter on the gravitational field of the string. 
We can also split the stress tensor into a static source, which includes the static part 
of the classical source and of the quantum stress tensor, and a time dependent part. 
The Riemann tensor is computed using some approximations in the stress tensor and 
we discuss specially its behavior at large times. The results are quite different if the 
string is wiggly or unperturbed. For wiggly strings there is a Newtonian potential at 
the classical level already and the quantum correction to such potential is too small 
to quantitatively modify wake formation behind a moving string. For unperturbed 
strings there is no Newtonian potential at the classical level but a Newtonian potential 
appears as a quantum correction and, thus, the string will exert a force on non 
relativistic matter. However this effect would be only significant when r is very 
small, say of the order of Tq, i.e. only at a microscopic scale. Perhaps back reaction 
might be significant to modify the dynamics of string crossing, however, we should 
keep in mind that when two strings approach at a distance of order vq the dynamics 
is dominated by the microscopic field dynamics and the effective macroscopic picture 
that we use breaks down. 
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2 Classical analysis 



In this section we compute the gravitational field created by the string as it is formed. 
We use the following classical stress tensor for the formation of a straight string which 
lies along the 2;-axis 



T/, = ^l5{x)5{y) diag(l, 0, 0, rd{t)), (2.1) 

where is the mass per unit length of the string, 6{t) the step function and r a 
parameter (0 < r < 1) used to modulate the tension, /ir, of the string. If r = 1 the 
string is an unperturbed cosmic string whereas if r< 1, this tensor gives a macroscopic 
description of a straigth but wiggly cosmic string. In this last case the effective 
equation of state for the string is p9|, |^ r/i^ = /iQ, where /io is the unperturbed mass 
per unit length {i.e. for a GUT string G/io~10~^). From numerical simulations one 
has in the matter era the typical value [Q /i~1.4/io {i.e. r~0.5). Note that when the 
string forms each segment may be approximated by an unperturbed cosmic string, 
a long segment becomes wiggly only after evolution of the string network by string 
intersection and by chopping off small loops. Thus back reaction is more important 
when an unperturbed string forms, since it takes the smallest time, however we keep 
T arbitrary in what follows since vacuum polarization will also exist once the string 
settles into a wiggly string and, in any way, keeping r arbitrary allows a simple 
identification ot the time dependent effects. 

Note also that in the classical stress tensor (|2.1|) we are using two simultaneous 
approximations. The first is the thin line approximation which assumes that the string 
has zero thickness, but, as we have emphasized earlier, a physical string has a radius 
To which gives a cut off radius inside which the approximation cannot be trusted. The 
second is the step approximation which assumes that the string is suddenly formed 
at time t = 0, the use of 6{t) instead of a smooth function which grows from zero to 
one in a certain time T requires a cut off in the momentum of the quantum modes 
of the order of 1/T. The time of formation T is bounded by the age of the universe 
when the string forms (for GUT strings T< lO^'^^sec). 

In the weak field approximation the metric in Cartesian coordinates can be written 

as 

gixy{x) = ViMu + V(a;), (2.2) 

where rjfj^u = diag(l, —1, —1, —1), | /i^j, |^ 1 and we have the gauge freedom, due to 
infinitessimal coordinate changes x'^ = x^ + C,'^{x), that 

for arbitrary ^^{x). Einstein's equations for the metric perturbation /i^j^ can be 
written in the harmonic gauge {h'^^—^^'n'^^h^a)^^ = as 

□ V = -IGTrGS^,, (2.4) 
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where S^^ = T^j, — \ rjfj^i.T'^^ and we still have the freedom within the harmonic gauge 
of choosing functions C,^{x) such that 0^f^[x) = 0. Note that a different analysis of 
this model which goes beyond the linear approximation has beeen given in ref. p8| . 



Since we start with a static prescribed metric before the string forms at t = we 
can find the solution of (|2.4|) for t > as a Cauchy problem. As it is well known, 
such a solution is given by Kirchhoff 's formula, 



h^,y{t,x) = / d^x' dt'h^,y{0, x') D{t, X — x') + hfj^i^{0, x') dtD{t, X — x') 

-WnGj^dt'J d^x' DR{t-t',x-x')S^^{t',x'), (2.5) 



where is the retarded Green function for a massless scalar field, i.e. DFt{x — x') = 
^S[{x-x'y]e{t-t') or -^S{t-t'- I x-x' I)/ I x-x' I, and D = Dr-Da is the 
Schwinger Green function {Da is the advanced Green function), which may be written 
as D{x-x') = ^5[{x-x'f] [e{t-t')-e(t' -t)]. The first integral in is over the 



hypersurface t' = and it is the solution of the homogenous equation which satisfies 
the boundary conditions given by the metric, h^y, and its first time derivative, dt'h^,,, 
at t' = 0. The second integral is the solution of the inhomogeneous equation which 
vanishes at t = {i.e. the boundary conditions are implemented with the solution 
of the homogeneous equation) and has support, see Dr, on the truncated past light 
cone starting at {t,x) and ending at the hypersurface t' = 0. 

The solution of ( |2.4| ) for the stress tensor ( |2.1| ) is somewhat simplified if we write 
( ^TTD as the sum of two stress tensors, both conserved, one of them, T^J,', static and 
the other, T^^J, time dependent, i.e. Tfj_y = Tj^^+Tj^^, with 

T(^) = ii5{x)5{y) diag(l, 0, 0, 0), T'^' = ^lT5{x)5{y) 9{t) diag(0, 0, 0, 1). (2.6) 

This leads to write h^u = h^^l + h^^l, where h^^l is the static metric corresponding to 
the source Tj^^ and is the time dependent part corresponding to the source Tj^y. 
The field equation (p.4|) for the static part h^^l is now simply 



V%1 = S'KG^x6{x)6{y) diag(l, 1, 1, 1) (2.7) 
which, using the cylindrical symmetry of the problem, has the solution 

h't^ = = = = 4Gf, ln(^ ] ^ a{r), (2 



where r^ = a;^ + ?/^ and R is an arbitrary constant with dimensions of length. 

The field equation for the time dependent metric perturbation /ijf^ is now given 

by 

nh^^l = -87rGfiTS{x)S{y) e{t) diag(-l, 1, 1, -1) (2.9) 



5 



whose solution as a Cauchy problem is given by (|2.5| ), but we note that since the 
boundary conditions for these components are obviously h'-^KO, x) = 0, dfh^^KO, x)=0 
we are left with the second integral of ( p.5| ) only. Using the second representation for 
the retarded Green function given above we have after a simple integration that the 
(unique) solution for is 

h'u' = h^^l = -K'l = -hf^ = 4G/irarccosh(^) e{t - r) ^ (3{t,r). (2.10) 

It is easy to check that the final metric perturbation h^y = h'~^l + h''^l satisfies the har- 
monic gauge condition. We can write the gravitational field of the string in cylindrical 
coordinates (t, r, 6*, z) as 

ds^ = {l + a + (3) d^-il-a- (3) dz^ - [1 - a + (3) {dr^ + r^dO^), (2.11) 

where a{r) and l3{t,r) are given by (|2.8| ) and ( p.lO| ) respectively. 

A few comments on the metric ( p. 11 ) are now in order. First, we see that the 



metric is continuous but its first derivatives are discontinuous along the light cone 
t = r. This is a consequence of the use of the step approximation: had we used a 
smooth time dependent function instead of 6{t) the derivatives would be smooth too. 
This means that we can only trust our results outside the spacetime region bounded 
hj t = r — T/2 and t = r + T/2. To this we should add that due to the thin line 
approximation the results are only reliable for r>ro, where is the string radius. 

Second, we also note that the metric perturbations diverge at r — ^ oo {t fixed) and 
at t — cxD (r fixed) but these divergences, as we will see from the Riemann tensor, are 
only gauge effects, they are a consequence of the use of the harmonic gauge. 

Third, one expects that when t ^ r {r fixed) the metric becomes that of a static 
cosmic string, unperturbed if r = 1 or wiggly otherwise. But this again cannot be 
seen directly in the harmonic gauge. The coordinate change which puts the above 
metric in a suitable form to see this is rather messy and, since it gives no further light 
on this issue, we shall not write it here. 

An important physical observable is the deficit angle of the two-surfaces /!: = const., 
2; = const, for the metric ( |2.11| ). Following ref. []T^, given a closed piecewise smooth 



curve on a two-surface which encloses a regular and simply connected region S, one 
defines the deficit angle associated with this curve as the angle A6 rotated by a vector 
after being parallel transported once around the curve. An application of the Gauss- 
Bonnet theorem shows that this angle is given by the surface integral of the Gaussian 
curvature K over 5*: A^^ = J^KdS. If the two-surface has circular symmetry the 
deficit angle is defined as the deficit angle associated with circles of radius r. It is easy 
to see that in a space-time described by a cylindrical symmetric metric the deficit 
angle of the two-surfaces t = const., z = const, reduces to 



Ae{t,r) = 27r 



1 -J—^ I 

V-9rr or 



[2.12) 
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hi -r hi A. Note that 



which in the hnear approximation becomes A6{t, r) = rr {h[. 
this definition differs from that used in ref. [|12|. For the metric ( |2.11|) we obtain 

t 



A9{t,r) = 47rG/i 



1 + r 



eit-r) 



2.13) 



It is clear that when t^r, r fixed, A6 ^AnG^ (1 + r), which is the deficit angle for 
a static string. 

The properties of the metric ( |2.11| ) are better deduced from the Riemann tensor 
which is gauge independent. Note also that some Riemann components give the tidal 
forces on surrounding test particles. In the linear approximation we have 



R 



i5',5p,d^ + d^Jld^dp - SPjp^d^ - d^d^dM hp,. 



[2.U) 



Following the separation of the metric perturbation into a static part and a time 
dependent part, h^^ = h'^^l + h'-^l, the Riemann components can also be separated as 
R^^lap and R^^lafj, the last is easily identified because it is proportional to r. For rj^O 
they are 



R 

Ra 

R 



tztz 



-R 



2GnT 



xyxy 



Vf 



xzyz 



R 



2Gfx 

Rtxty 

2Gi2 



cos 20 
2G/i 



-rVf 
sin 20 



(62-1)3/2 

b 



e{b-i] 



1 - 



cos 20 



cos^ 9 



T 

-Vf 
2 ■' 



(62 

6(262-3) 



0(6-1) 



(62-1)3/2 



0(6-1) 



txtx 



R 



tzxz 



R 



cos20-rP/ 
2G/ir 



tyyx 



COS0P/ 



COS 20 + 



sin^ 



(62-1)3/2 



(62-1 
0(6-1) 



0(6-1) 



(2.15) 



where we have introduced the new variable 6 = t/r, instead of t, and where Vf 
denotes the Hadamard finite part, which gives well defined expressions, in the sense of 
distributions, on the light cone 6=1 (see Appendix A). Ryzyz, Rtyty, Rtzyz and Rtxxy 
can be obtained from Rxzxz, Rtxtx, Rtzxz and Rtyyx respectively by interchanging 
COS0 and sin0. 

It is now clear that when r = 1 (unperturbed string) the Riemann tensor vanishes 
when t ^ oo (r fixed), i.e. limt^oo Rfiua/3 = 0, and the spacetime becomes flat. But 
if r 7^ 1 there are tidal forces among test particles wich correspond to a Newtonian 
like potential /i^f oclnr. Note that a Riemann component such as Rtztz, which gives 
relative accelerations among particles along the direction of the string, is a transient 
term, i.e. Rtztz = when t <r (in regions not yet affected by the formation of the 
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string) and Rtztz — 2Gfir r^"^ b^^ (l + | + O {b^^)j when t>r, which approaches 
zero very quickly. 

Our next task is to obtain the quantum correction to this curvature tensor due to 
the quantum fluctuations of matter fields. 



3 The stress tensor for matter fields 



Quantum fluctuations of matter fields interacting with the gravitational field of a 
cosmic string give a non null vacuum expectation value for the stress tensor of these 
matter fields, even if they are conformally coupled. For a free massless conformally 
coupled scalar field in a flat spacetime background with arbitrary linear gravitational 
perturbations ( |2.2| ) this stress tensor has been computed to one loop order by several 
|2^, |23, E5|, E6[ and it is given, to first order in hu^i^x), by 



authors 



{T>^^(x)) = -— B'^^x) + 3ah J d'y H{x-y, fi) A^-'iy), 



(3.1) 



where a = {288QTi'^ 



B 



flU. 



x) = 2r/''^nG/ - 2G/''^'^ 



2 2 
-G "•^'^ + -v^^aG " 



A^^'ix) 

G^'^ is the Einstein tensor for the metric to first order in h 
propagator defined by 

1 /■ d^p 



(3.2) 

^1., and H{x-y, ft) is a 



H{x-y,iJ,) 



-ip(x-y) 



In 



It 



+ 2Tiie{p^)e{-p^ 



(3.3) 



where Jx is an arbitrary renormalization scale. Notice that the second term in (|3.1|) is 
traceless and that it is the term proportional to B^^^{x) which gives the trace anomaly 
in this case, 

ahUR, (3.4) 



(T';(a;)) = -— E'^, 



x] 



where R is the scalar curvature to first order in /i^^,. Using Vfj.G^" = 0, it is easily 



seen from ( p.2|) that the two terms in (|3.1j ) are independently conserved. 

We can now proceed to the computation of {T^^{x)) to first order in h^j^. In 
section 3.1 we will calculate it outside the string, that is, for r^Q. In section 3.2 we 
will see how a generalization to include r = can be made. In section 3.3 we give 
the approximation to the stress tensor which will be finally used to compute the back 
reaction on the string metric. 
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3.1 The quantum stress tensor outside the string 



From ( |3.1|) and ( p.2|) we see that all the dependence on the metric perturbations is 
in the Einstein tensor 0^,^. Now since the source of such gravitational perturba- 
tions is the classical string source T^^'^ given in we can use Einstein's equations 
G'^'^ = —SttGT^'^'^ in ( |3.2| ). Thus it is worth to remark that we do not need to know 
the explicit gravitational field created by the forming string, it suffices to know the 
explicit form of the classical stress tensor which produces that field. After the above 
substitution, since T^'^'^ is proportional to 6{x)6{y), both B^^{x) and A^^{x) have 
support on t he string core. If we are interested in computing {T^^{x)) outside the 
string, i.e. for tt^O, it is clear from ( |3.1| ) that the only contribution comes from the 
second term in ( |3.1D and that terms in H{x — y,fi) with support on x = y will not 
contribute. 

Let us give a more suitable representation of H{x — y,Ji) for x^y. For this it is 
useful to define the four-vector 



-ip{x-y) 



d 



In 



'-(p^ + ie) 



(27r)4 5p« 
This vector can, on the one hand, be simply computed as 



27ri0(/)^(-j9°) 



(3.5) 



Fa{x-y) 



d 



(9x° J (27r)^ 



-ip{x-y) 



1 



+ ie 



2ni6{py{-p° 



d 



dx' 



-Dnix-y) = -ix-y)^6'\ix-yf ^^(x° 



TT 



y 



(3.6) 



where S' means derivative with respect to the argument of 6 and we have used a well 
known integral representation of the retarded Green function Dji and the form of Dji 
given in section 2. 

On the other hand, integrating (|3.5|) by parts, we can also write Fa{x — y) as. 



Now comparing (|3.6|) and 



Fa{x-y) = i{x -y)aH{x-y,ij). (3.7) 
we get the following useful representation for H, when 



1 



H{x-y,jj,) = -6' (x-yf ^(x° 
Thus using ( |3.8| ) we have outside the string 

3ah 



y 



(T'^ix)) 



TC 



with 



d'yS' {x-yf e{x'-y')A'^''{y), 



(3.8) 

(3.9) 
(3.10) 
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where we have defined 



/ per 



If we now write, see ref. 



5' 



{x-yf 



{x-yf + X 



in (|3.9| ), change variables and use 



= hm 6' 
we have exphcitly 



: hm —5 

A^o- dA 



(3.11) 



d_d_ d_ 

p^'dx^dx^ A™- dX 



d^x' 5{x'^ + X)e{t') 



x5(x - x')5{y - y') {5151 - r9it - t') 5^51) (3.12) 



and the computation is now straightforward. It is clear from (|3.12| ) that we have 
two types of terms in the stress tensor, (T^^) 



. 7^(1) 



(2) 



the static terms, T, 



(1) 



''A^lu^ ^^^^ otduiu uciiiio, -LA^lIy^ 

which depend on the static part of T^^'^ and the time dependent terms, T^^^i,, which 
depend on the time dependent part, that is, the terms which are proportional to 
T6{t — t'). The sub index a stands for the part of the stress tensor related to the 
tensor A'^'^ in ( p.l|) . The static terms are easily computed, the / dx' dy' integration is 
trivial and we use 5{x''^ + A) to perform the / dt' integration, following ref. [^|. These 
terms are all proportional to The time dependent terms, although sligthly more 
complicated, can be computed without difficulty. The final result for r 7^ in the 
Cartesian coordinate basis is: 



2a 



[4-r/i8(6)], {T:) 



(T,'^) = — — r^7(fe). 



4(T 



12 



TO = - 3r [2 + r/i2(6)]-4 



XX 



[2 + r/i5(6)]|, 
(3.13) 



where a = SaGfih, the index a = 1,2 refers in to x^ ■ 
coordinates), and 



-X, X 



-- y (the string transversal 



9{b) = Vf 



1 



(52 _ 1)5/2 



e{b-i) 



Mb) ^ Vf 



(P2 _ 1)5/2 



(26^ -562 + C/4) 0(6-1) 

(3.1' 



As above Vf denotes the Hadamard finite part, which gives well defined distributions 
on the light cone 6 = 1 (see Appendix A). It is easy to identify directly from (|3.13|) the 
static and time dependent parts of the stress tensor, i.e. T^l^i, and T^l^^ respectively. 
Note that T^l^^ corresponds to the vacuum expectation value of the stress tensor 
due to a dust rod along the 2;-axis with mass per unit length /x. It can be checked. 



computing the derivatives of the distributions ( |3.14|) using the methods of Appendix 
A, that the (T^i/) just derived is conserved. 

The most salient feature of (T^v) is that it quickly settles down to the final static 
values. Note that when t grows keeping r fixed, b = t/r grows as t and we can expand 
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g{b) and fc{b) in terms of r^: ^(6) -6^5 + | 6"2 + (6"^)), /c(&) ~2 (1 + (6^^)). 
This means that fc'{b) differs from the static value, 2, by terms of order and that 
g{b) goes hke b^^. Thus (T^'^(x)) is effectively time dependent only when t^r, it 
reaches the static values very quickly. 

The final static values of (T^jy), i.e. when t —* oo, are read off from ( |3.13|) . In the 
polar coordinate basis {dt, dr,do,dz) they are 



t— >oo 



diag(4 - 2r, 1 + r, -3 - 3r, -2 + 4r). 



(3.15) 



In particular when r : 
ofrefs. EOlBTl. 



1, i.e. the string is not wiggly, we get the well known results 



3.2 The quantum stress tensor including r = 

Since the quantum stress tensor is a source in the semiclassical Einstein's equations, 
we need to know {T^u) in all the space-time. One way to do this is to try to generahze 
the previous calculation to include r = 0. Such a generalization should have the form 



(T'^^(x)) 



Trix)+T, 



X) 



(3.16) 



where Tg"^ 



-{ah/ 6) B'^^ and T^^^'^'^ and T/^^^^ are the traceless tensors corresponding 
to the static and time dependent parts respectively in the second term of ( |3.1| ). These 
should be some well defined distributions which reduce to the expressions (|3.13| ) for 
rT^O. Using Einstein's equations G^^ = in ( p. 21 ), Tg"^ can be expressed in 

terms of the classical stress tensor so its computation is straightforward. One finds 
the following non-null components in the Cartesian coordinate basis: 



-a{l + Te{t))V' {6{x)6{y)) 



Tst = -oT5{t)da{5{x)5{y)) 



Bz 



rp a 

'''Bb 



TT 

-— a 
3 

TT 



2 crj^ab 
-fl 



1 + T9{t)) iS{x)6iy)) -r5'{t) 5{x)5{y) 
1 + Te{t)) {5{x)5{y)) -r5'{t) 5{x)5{y) 
+ Te{t))dad,{5{x)5{y))}. 



(3.17) 



It is easy to see that the tensor Tg just found is conserved. From the expressions 

refers 



( p.l7| ) it is clear that we can write j'^'^ = j'^^^t^'^ ^j'^^^t^" where, as always, T^^ 
to the static terms and T^^^^^ to the time dependent terms. 

To generalize the calculation of the second term in ( p.l| ) to include r = 0, we need 
to extend the representation (|3l8| ) of the propagator H{x—y, Ji) to all values of {x—y)". 
Such an extension, which is derived in Appendix A, is 



H{x—y,fi)= lim^j— 5' [x — yY 6{x 







1/°) 9 i\x-y\ -e) + [In/ie + 7-I] 5\x - y) 

(3.18) 
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where 7 is Euler's constant. For the static terms T^''^^ we find 

fj^* = -47raV2/(x,y), fj^^ = 27raV2/(a;,y), 

fi^'t = 27r(T (5,5V^ - dad,) lix.y), (3.19) 



where I{x,y) = J d^x' H{—x',fi)6{x + x')5{y + y'). This integral is computed in 
Appendix A with the result 



£m |i- l^^(r - e) + (in f + 7) 5(a:)<5(l/)} • (3.20) 



Using the result given in ref. V P/ (r ) =4P/ (r ) — 27rV (5(x)(5(y)) one can 
substitute in these expressions 

V'/ = ^ ^/ (;^) + ^ + ^ - 1) V'(5(x)5(y)) . (3.21) 

These results for T^^^^'^ agree with the static components of (|3.13|) when r^O. 

For the time dependent terms T^^^'^ the calculation is considerably more com- 
plicated. Thus, instead of computing these exactly, we will introduce the following 
approximat ions . 

3.3 The approximated quantum stress tensor Tj^^^" 

We introduce here some approximations in T^^^^ . First, instead of working out the 
exact expresions for these terms to include r = 0, we will introduce a cut-off radius 
To (which can be viewed as the physical string radius) and assume that the tensor 
j^miJ'i' (jgriygfj in ( I3.13D is only valid for r > tq. We will make some assumptions for 



the values of this tensor at r < tq. At the end of the calculations we will take the 
limit To ^ 0. Beside this, in the terms of ( p.l3|) which correspond to T^^^'^ we will 
substitute the distributions fcip) and g{h) by 2 9{t—r) and zero respectively. That is, 
we make a sudden approximation assuming that the values of T^'^^^ change suddenly 
at the light cone t = r. This seems justified in view of the fact that such terms, as 
we have seen before, settle quickly to the final static values. But now, to ensure the 
conservation of T^^^^^ on the light cone t = r we need to add terms proportional to 
6{t—r), which have the same singular behavior on the light cone as the stress tensor 
( p.l3| ). Imposing also that T^^'^^ = we obtain the following approximated stress 
tensor for r > ro : 



2ar 



/ 2 



-4 ) 



( 1 
1 

V 



-3 / 



(3.22) 
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in polar coordinates. Now we have to make an assumption for the values of T^fl^, at 
r < To- We take a stress tensor with terms proportional to 6(t — r) and r6{t — r) as 

)nstant for t> r. Im 
0, we find for r<ro: 



before, but which is constant for t> r. Imposing that it is conserved, '^JT^^^^ = 



and traceless. 



-'-A 



rp (2) 



2crr 



Bit- 



i 



\ 



\ 



+ r5{t- 



( 



V 



Note that when we take the divergence of the complete Tj^^^^ we have the step 



-1 / 



. (3.23) 



functions 9{rQ- 
derived too. 



and 9{r — rQ) multiplying the above expressions which have to be 



It can be seen also that the values of ^ 



correspond to the values at r : 



the classical stress tensor for the string given in refs. [|T5| [16 



for r < ro and t>r obtained in ( p.23| ) 
of the quantum stress tensor for a model based on 

In such a model the 



string has a finite radius tq inside which the classical stress tensor is assumed to be 
constant; this tensor is given by 



^^'^^ =^^('^o-r)diag(l,0,0,r). 



(3.24) 



where we have introduced the r parameter in order to allow for the possibility of a 
straight wiggly string. The string energy density e is related to its energy per unit 
length /i (see refs. ll^, |l6l) by 4Gfi = 1 — cos (^roVSirGej. At first order in Gn this 
relation gives 7iGerQ = Gfi + 0{G'^fi'^). Using the previous expressions (p.l|)-(|3^), we 
can calculate the stress tensor (T^u) to first order in Gfi for a free massless conformally 
coupled scalar field with the classical source ( |3.24D . The values of such tensor at r = 
agree with these of the 6{t — r) terms in ( |3.23|) and this gives further justification to 
the approximation taken. 



4 The back-reaction metric 

In this section we compute the correction to the gravitational field due to the vacuum 
polarization of matter fields at one-loop order given in the previous section. For this 
we solve the semiclassical correction to Einstein's equations, 

G^'^(x) = -SvrG [T/''(x) + (r^"(x))], (4.1) 

in the linear approximation. Within this approximation the metric tensor can be 
written as 

gf^u = V^iu + hf,y + hf,y, (4.2) 

where h^j^y is the metric perturbation due to the classical string stress tensor TJ^'^ and 
h^y is the quantum correction to the metric which is of order Ti. The classical part. 
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h^y, has already been computed in section 2, and we note here that {T^^{x)) depends 
on that part only, since including h^^j in (T^^^x)) would lead to terms of order fi^ 
which we neglect. Thus eq. ( [4.1| ) leads in the harmonic gauge {h'^'^—^ri'^'^h"^),,^= 
to the equation 

□V = -IGnGS^,, (4.3) 

where S^i, = (T^u) ~ ^Vf^uiT'^a)- formal solution of this equation is given as an 
initial value problem by (|2.5| ), where one substitutes h^u by h^i, and S^i, by 5"^,^. Now 
from the decomposition (|3.16|) for the stress tensor and the separation of Tb^u into a 
static and a time dependent parts, '^b^^.u^ is clear that a similar separation 

can be performed in S^^, i.e. 

S^v = T}^Iu + "S^/li/ + '^A^u + ^bIu^ (4-4) 



where we use that Ta/^i/ is traceless. Then ( [4.3| ) leads to likewise separate quantum 



perturbations: hfj,u = h^j^^^ + h^lti, + h^lu + ^i^/ii/ = ^''/ll + ^JS- The main advantage of 
this separation is that the initial conditions at the surface t = for h^^l are simply 



"'fJ.L' 



= 0, (4.5) 



consequently, the explicit form of h^J^^^(t,x) is given by 

hi'Ut,^) = -ISvrG /* dt' f d'x' Dn{x - x') f^^^H', x'). (4.6) 

Similarly hj^^^,^ is given by this expression changing T^^^^^ by S^^^^. 

We are only interested in finding these semiclassical perturbations h^^, outside the 
string, that is, for r^Q. 



4.1 The static part 

For the static correction h^^l eq. ( |4.3| ) reduces to 

W'h}^ = IQt^G S% (4.7) 

where S'^J,^ = T^^^^ + 5^^^,^. It is very simple to solve these equations for r 7^ mak- 
ing use of the cylindrical symmetry of the problem. Notice that 5^^' has support 



on r = 0, so from ( 3.13|) , we have in the polar coordinate basis S^^^^ = T^^j^^ 



—2(jr~^ diag(4, 1, —3, —2) for r 7^0. It is now easy to see that the following functions 
(which are defined in terms of the metric perturbations), = h''^^l, = h^'^'^l, 
h^^^ = h^'^^l + h^^^Q and h^^'' = h^^^^ — h^^^g, depend on r only and satisfy the equations, 
for ry^O, 

GAnaG 
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where 71 =74 = —2 and 72 = 73 = 1- The harmonic gauge condition in terms of these 
functions can be written as (/i^^' + /12 ' — ^4^'),r = 2r~^hll\ The general solution of 
these differential equations can be expressed as a linear combination of terms such as 
r^, Inr, and r~^lnr. Finally, after a gauge transformation, within the harmonic 
gauge, the solution can be written as 



(1) 



-327TaG — + aG A In 



r 



= leTraG— , 



aGA In 



1 



r 

R'2 



32naG— ln( - 
R 



(4.9) 



where the constants of integration Ri, R2 and R have dimensions of length and A has 
dimensions of (length)"^. Notice that there is a gauge freedom within the harmonic 
gauge to set the constants Ri, R2 and i? to a fixed value, ie. a c hange in the values 
of these constants is equivalent to an harmonic gauge transformation. The values of 
these constants in the solution for r 7^ can be determined if we solve the equations 
( [4. 71 ) including r = 0. However, without working out the explicit solution, we can use 
dimensional arguments to see that A must be zero. Note that the only dimensional 
constant parameter which can be used to make h^^l/aG is /i, but since the dependence 
in this parameter must be logarithmic we are forced to take A = 0, on the other hand 
R will be proportional to In fact, we can set R=l/fi after an harmonic gauge 
transformation for r^O. 

By doing a little more work, we can also arrive at the above solutions solving 
explicitly equations ( [4.7] ) with the inclusion of r = 0. In fact, from ( |3.17D and ( ^.19| ) 
it is easy to see that the equations for hi \ /ig' and /I'^'i^ + Zi*^'^ are 



AI--6{x)6{y) 
o 

21 + ^6{x)6{y) 



0, 



32n^(xG 



I+-S{x)S{y) 



0, 



(4.10) 



where / is given by (pOD. These give the solutions h[''> = -IGn^aG (4 / - | 6{x)5{y) 

:327rVG' (/ + l6{x)6{y)), which for 



levrVG (2 / + I S{x)6{y)) and /I'^'l^ + Zi'^'^ : 



r 7^ reduce to 



with A = 0. If we write the remaining two equations for 



/I'^'i^ — /i*^-'^ and h^'''y, using that S{x)5{y) = 2^ V^lnr, we see that for r 7^ one has 
h''^^ = 327raGr~^ ln(/t/2r), where k is some numerical factor. So we finally have for 
r^O 



(1) 



-32'KaG^, 
IGiraG — , 



IGvraG — , 
32'naG — In Jir. 



(4.11) 



This solution can also be obtained by the procedure which will use next to obtain 
the time dependent metric perturbations. Namely, we introduce a cut-off radius 
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ro and take the limit tq ^ at the end of the calculation. For this we use an 
approximated stress tensor for r<ro as in section 3.3. 

Let us now consider the contribution A6"^^ to the deficit angle of the two-surfaces 
t = const., z = const, due to these static corrections to the metric. Substituting in 
( p.l2| ) the values of h^^l given by ( [4.11| ) we find 

A^~« = (4.12) 

We can now calculate the static quantum correction to the Riemann tensor. Using 
the expressions derived in Appendix B we find the following non-vanishing compo- 
nents for r 7^0 

IGnaG , „ , 327ro"G „ "ilixoG 

R'^L. = (2cos2^ + 1), = sm2e, R%,y = — 

R[L = (2 COS 2^ + 1), R'll,y = sin 2^, (4.13) 

and also R^^lyz and R%ty^ which are obtained from -R^^^xz R^tltx interchanging 
cos^^ by sin^^. 



4.2 The time dependent part 



Let us first consider the contribution to the quantum perturbations to the metric 
coming from T^^l^^ given by ( p. 22 ) and (p. 23 ). This solution is given by 



(2) 



Defining, 

as we have done before, h^^^=h';^^\, h^^ = h'^^l, hj^^ = h'^^l.+h^^l, = ^A^r~^A'e 
hj^^ = h^^^, this leads for (t— r)>2ro and r>ro to the solution 



(2) 



32aTG r 



ddgk{d) 



Jo 



Xo 



dx X 



+ 



-1 



2 \ 6— cos 6 



Gk + CkX 



dx 



d_ 
db 



d 

Dk + DkX — 
oh 



fib;x,9), (4.14) 



where b = t/r, xo = ro/r and 
Gk = —2 6k^ + 4:6k2 — 2 6, 



fc = Ofcl + 4 Sk2 - ^ <Jk3, 

-Dfc = 2 4i - 4 5k2 + 2 5fe3 - 4 5k4, 
1, A; = 1,2,3 

cos 2^, fc = 4 
cos 6*, k = 5, 



9kiO) E 
f{b-x, 9) = In 



Gk = Ski + (^^2 — 2 Sk3 + Sk5, 

Dk = 4i - 3 6k2 + 2 5fe3 - 4 6k4 + 45, 



b-x + d2x{cose - 6) + 62 _ 1 



ln(x^ - 2xcos6' + 1). (4.15) 

2 
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We now expand the solution ( [4 .141) in terms oi XQ = ro/r neglecting terms of order 
Xq wich will vanish in the limit ro— *0. We find 



h 



(2) 



327rarG' 



(4.16) 



where 



= = 



2 (62-1)3/2' 



$" = 0, 



2 (62 _ 1)3/2' 



1 

1 + - 



(62 - 1) 

6 

2 (62-1)3/2^ 



+ 



(62 - 1)3/2 



= :t -1 + 



(62 
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(266 _ 2^2 _ 3) 1 



(62 - 1)3/2 

6 



(62 - 1) 



4 (62 - 1)3/2 



- + arccosh 6 — In 4(6^ — 1) 
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— + arccosh 6 — In 4(6^ — 1) 



(26^ - 1) (46^ - 11) + arccosh 6 



1 



2 (62 - 1)3/2 



(26^ - 3) 



3 In 26 - 2 arccosh 6 - In 4(6^ - 1) 



^1 = 



5 1 
■4 62 
5 



3^2(1, l,^;2,4;i 



(62 - 1) (62 - 1)3/2 



+ 62 + arccosh 6 - In 4(6^ - 1) 



(4.17) 



It is easy to see that we can change in ( [4.16| ) lna;o by ln(i?/r), where R is an arbitrary 
constant with dimensions of length. For this we make the following harmonic gauge 
transformation 



it = —^'kotG — , In 
£x = SnarG — ^'^ In 



0, 



We can also set -R=l//i and, after this, we can take the limit tq-^O in hj^ 



(2) 



Let us now consider the terms /ie^^,. Substituting T^^l^, by 5|b^^ in 
defining h^]^ as before we find for r 7^ 



and 



L (2) 



327r(rTG 1 

J.2 Q 



vf[<!>mo{b-i) 



(4.19) 
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where the functions $^(6) are defined in (|4.17|) . It is easy to see that these terms can 
be ehminated with a gauge transformation of type ([4.18| ). Thus we finally have for 
r 7^ and 6 > 1 

ht' = ^ mib) - In/ir $^,'(6)] . (4.20) 

We can now compute the time dependent quantum contribution A^^*^' to the 
deficit angle of the two-surfaces t = const., z = const. Substituting in ( p.l2| ) the 
previous values of h'^^l we get 



ib' - 1) 



(46^ + 136^ - 2) + 



36 



(62 - 1)5/2 



- (86^ + 86^ - 15) 
4 



5 1 



+3 (ln26-arccosh6) - - — 3F2 ( 1, 1, -; 2, 4, 



4 62 



7 



62 



(4.21) 



Expanding this expression in t ^ with r fixed one finds 



(2) 



(4.22) 



that is, A6"^' reaches its static value 327[^o'tG r""^ very quickly. Note that this con- 
tribution due to the time dependent terms is exactly the same as the contribution 
A^^'^^ of ( 4.12 ) when the string is unperturbed (r = l). 



Finally we can calculate the time dependent part of the Riemann tensor 



The substitution of the terms hk of Appendix B by the terms h'j^^ of ( |4.20[ ) gives the 
following result for r 7^ and 6 > 1 



2ATTaTG 



tztz 



R 



(2) 

xzxz 



xzyz 



txtx 



6(2) 
^txty 



r4 (62 - 1)7/2 



(26^ + 3), 



A-narG 



inarG 



A-narG 



+ 



36 



(62 - 1)7/^ 



(26^ + 3) + COS 2^ 16 + 



156 



(62 - 1)7/2 



sin 29 16 + 



36 



15 6 



(62-1)7/2^ 



(62 - 1)7/2 



(26^ + 3) + cos 26 



+ 



15 6^ 



(62 



(26^ + 3) 



15 6 



(62 - 1)7/2 



i (86^ + 86^ - 11) + 3 (ln26 - arccosh6) 



A-narG 



sin 2^ <^ 8 + 



15 6^ 



(62 



5 1 
~4 62 

(26^ + 3) - 



3i"2(l,l 



7 



5-^)0''' 



2,4; 



62 



156 

(62 - 1)7/2 



^6^ + 86^ - 11) 
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+3 (In 2b - arccosh 6) - ^ 1 3F2 (^1, 1, ^; 2, 4; ^ 



RfLz = ^^^^^^ cose 



1 



5(2) _ l^TTOrrG 
J^txxy — 1 bin (7 



xyxy 



AnarG 



{p2 _ 1)7/2 
\ (62 - 1)3 



(46^ + 1), 

(86^ + 1452 + 3) 



(4fe2 + 1) 

(62 - 1)7/2 



4 



(86^ + 86^ - 15) 



+3 (In 26 - arccosh 6) - ^ i 3F2 (^1, 1, ^; 2, 4; ^ 



1 



(62 - 1; 



(166^ + 246^ + 396^ - 4) 



36 



+ (6^3T)V^ + 2) 



^ (86^ + 86^ - 15) + 3 (In 26 - arccosh 6) 



5 1 
'4 62 



— ^3i^2(l,l,-;2,4;^ 



7 

2' 



1 

62 



(4.23) 



where 3F2 denotes a generahzed hypergeometric function. As before R^ylyzi R%tyi 
R'tzyz and B}tyyx are obtained interchanging cos^ by sin^ in R^^^xz^ R'tltx-, R^tLz and 
R^'txxy respectively. Notice that R^^i,a/3 does not depend on the arbitrary parameter fi. 
This is due to the fact that, as we have seen before, the value of /2 can be changed by a 



gauge transformation. The static limit of these components, i.e. limi? 



(2) 



= R' 



IS 



T3L 

xzxz 



R 



txtx 



327rcrrG 



;2cos2^+l), R 



L 

xzyz 



(2cos2^^ + l), R\,,y 



QATxarG 
327rarG 



sin 26, R 
sin 29, 



xyxy 



327rarG' 



(4.24) 



and R^yzyz an d R\y ty are found interchanging cos 6 by sin 6 in R^^^xz and R^^t^ respec- 
tively. From (|]2|) we see that these final static values are quickly reached since the 
corrections for large times go at least as 6~^ = r^/t^. The final semiclassical Riemann 
components are obtained by adding to the classical values ( |2.15D the back-reaction 
corrections ( |4.13D and ( |4.23| ). 



5 Conclusions 

In this paper we have derived the back reaction, due to quantum fluctuations of matter 
fields, on the gravitational field of a cosmic string during and after its formation. As 
matter fields we have just considered a massless conformally coupled scalar field but 
the results are easily extrapolated to of such fields. As a model of cosmic string 
formation we take an initial thin straight rod whose tension grows suddenly (step 
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approximation) from zero to a maximum value which corresponds to the mass per 
unit length if the string is unperturbed, or to a smaller tension if the string is wiggly. 

Within the linear approximation of Einstein's equations we have first computed 
the metric perturbations and the corresponding curvature tensor for both unperturbed 
and wiggly strings. If the final string is unperturbed the Newtonian potential vanishes 
and the spacetime becomes flat with a deficit angle, but if the string is wiggly there 
is a Newtonian force per unit mass which goes like r)r"^. This force may play 

an important role in the formation of wakes behind long strings |Q. 

We have then computed the vacuum expectation value of the stress tensor of 
matter coupled linearly to the string gravitational field. We have seen that after 
formation the stress tensor settles quickly to that of the final unperturbed or wiggly 
static string. The stress tensor has an energy density which goes like NhG^r~^ (as 
typical of Casimir type energies), where we have assumed matter fields which 
correspond to massless fields or with masses less than ntH, where tuh is the Higgs 
mass responsible for symmetry breaking. 

With such a stress tensor as a source we have computed the perturbation induced 
on the gravitational field. We have also computed the Riemann tensor components, 
which give essentially the tidal forces, corresponding to such quantum corrections. 
The static part of this tensor correspond to a Newtonian force per unit mass which 
goes like NhG'^fir~^. The time dependent part of the semiclassical perturbations to 
the metric and the Riemann tensor, on the other hand, are quite complicated. But, 
for a given radius, after a short time t>r the curvature tensor becomes static. Unlike 
for the classical part the quantum correction to the gravitational field does not differ 
substantially if the final string is unperturbed or wiggly. 

Let us now discuss the importance of the quantum corrections in the case of 
a wiggly string and the unperturbed string after their formation. In the case of 
a wiggly string, i.e. r < 1 (typically r ~ 0.5) [||, the ratio between the quantum, 
Fq and classical, Fc forces on surrounding nonrelativistic particles can be estimated 
in the static limit as Fg/Fc ^ GAira NhG {2 — t)/ [{1 — t) r"^], i.e., if Ip is the Planck 
length, Fg/F^ ~ IQ-^N (Ip/ry. This means that, unless is unreasonably large, 
the quantum effects are always negligible if r > 10/p. Note that for a cosmic string 
the smallest value that r can take is tq ~ h/mn and, since /x ~ mjj/h, we have 
Fg/Fc{r) < Fg/Fc (ro) ~ lO-^A^G/i which for GUT strings is of order lO^^ - lO^^. 

For an unperturbed string r = l there is no classical force Fc. Thus the quantum 
correction will be responsible for a Newtonian force on the classical matter surround- 
ing the string. However, that force decrases like which means that it becomes 
negligible very quickly at macroscopic distances from the string. In this case it is 
better to consider the deficit angle which is an important physical observable. The 
deficit angle can be written as a sum of a classical term plus a term of quantum 
origin + A^, using the notation of the last section. As we have seen, for the 
classical part AO^Gfi whereas for the quantum part A6' ~ 10~^A^/iG^/ir~^, so that 
Ae/A9 ~ IQ-^N {Ip/ry which is of the same order as the ratio between the quan- 
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turn and classical forces given before. The ratio of the transverse velocities on the 
string wakes due to the quantum effect and the classical deficit angle are also of this 
order. These results hint that whereas one should not expect any quantum effect at 
macroscopic distances on surrounding matter, quantum effects might start to become 
important near the string at microscopic distances. In particular, when two strings 
cross, there might be a correction, perhaps in the form of a Casimir like force, due 
to quantum effects. However in this case our classical picture breaks down and one 
must consider the dynamics of the Higgs fields themselves at microscopical level. 
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Appendix A 



A.l The Hadamard finite part 

We have introduced in this paper some singular distributions denoted by the symbol 
Vf. They are generated by the Hadamard finite part of a divergent integral. See refs. 
|]35| , |36| for more details on these distributions. The idea is the following, suppose 
that we have a function h{x) which is Lebesgue integrable on all the intervals (a+e, b), 
with e>0, but which is not integrable on (a, b) (a, b are some real constants). Let us 
consider the integral 

rb 

J(e) = / dxh{x), (A.l) 

Ja+e 

which diverges in the limit e ^ 0^, and assume that it is possible to separate it in 
two parts 

J(e) = /(e) + F(e), (A.2) 

where J(e) is a finite linear combination of negative powers of e and positive powers of 
In e and F{e) has a finite well defined limit as e — > 0^ (in the sense that is independent 
of the way by which is obtained). Then the Hadamard finite part of the divergent 
integral dx h{x) is defined by 

J^^r dx h{x) = lim F(e) = lim | dx h{x) - /(e) I , (A.3) 

J a 6^0+ e-+0+ [Ja+e J 

that is, we throw away the divergent terms in ( [A.2|) and then take the limit e ^ 0^. 
The definition is readily extended to functions h{x) defined on all the real axis which 
have singular points. All we have to do is to decompose the integral over iR as a sum 
of integrals of the kind considered above. Then one can define distributions denoted 
by Vf [h{x)] in the following way 

dx Vf [h{x)] (p{x) = J-p dxh{x)(p{x), (A.4) 



where (p{x) is an arbitrary tempered test function. The definition can be easily 
generalized to the case of several variables whenever the divergent integrals can be 
reduced to one-dimensional divergent integrals. 

As an example we calculate the derivative of the distribution (6^ — 1)^^/^ ^^(6—1) in 
detail. Notice that it is a distribution in a one- dimensional space (and also in a four- 
dimensional space if we set b = t/r) since the integral db (6^—1)^^/^ 9(b—l) Lp{b) = 

db (6^ — 1)"-'^/^ Lp(b) is convergent. To start with, it is easy to show that 

^ [{b^ - 1)-^/' e{b - 1)] = Vf [-b {b^ - 1)-^/" e{b - 1)] . (A.5) 
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For this, let us consider the following integral 

d 



dh 



l+e do 



roo 

where we have integrated by parts in the last step. This shows that 



dbVf \-b{b^-l)-^/H{b~l)] i^{b) = J^pi db 



db 



db (6=^-1)-^/" ^\b) 
d 



db 



{b-'-iy'lHib-l) ^{b), (A.7) 



where the definition of the derivative of a distribution has been used, and this proves 



( |A.5| ). An explicit expression for the distribution Vf \-b {b'^ -l)-^^"^ e{b-l) 
derived as a distributional limit. In fact, from ( |A.6[) we have 



can be 



dbVf -b{b^-l)-^/^e{b-l) ^{b) 



lim 

e^0+ \Jl+ 



lim 

e^0+ I Jl 



1 



lim 



db [-6(62-i)-3/2 

' db [-6(62_l)-3/2 

-6(62-l)-3/2^(6_l-e) +'s{b-l) 



b=l- 



db 



which gives the sought expression 

d_ 

db 



- i)-V2^(5_ 1) =p/ -b{b''-i)-''/''e{b-i) 



-3/2 



lim 

e-*0+ 



1 



-b {b^ - 1)-^/' e{b- 1 -e) + ^ (5(6 - 1) 



(A.9) 



where this distributional limit has to be understood in the sense of ( |A.8D . 
In a similar way it is easy to show that 



db^ 



(,2 _ iyl/2 Q^^ _ _ ^(26^ + 1) (62 _ l)-5/2 _ 

(26' + l)(62-l)-5/2 0(6_l-e) 



lim 



1 1 1 

+ 



,(2e)3/2 8 V2i 



5{b-l) + ^^J\b-l) 



(A.IO) 
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Other Hadamard finite part distributions appearing in this paper are given by the 
following distributional limits: 



— lim 



1 



Vf \np{b) (62-i)-3/2 ^(6-1)1 = Kp{b)Vf [(62-l)-3/2^(6-i) 



= lim 

6^0+ 



g{b)^vf\{b'-i)-'/'e{b-i) 



■ lim 

e-*0+ 



Kpib) {b^-l)-^/^eib-l-e)-K^5ib-l) 



fcib) =Vf [b (62-i)-5/2 ^2b^-5b''+C/4) 9{b-l] 

= lim|6(62-i)-V2(26^-562+C/4)e(6-l-e) 



+ 



C\ 1 1 fC\ 1 

1-— I T:r^ + - 5 + - 



12 J (2e)3/2 ' 8 V 4; 



2 V 127 V2i 



where /«7^0 is a constant and p{b) is a polynomial in b with p(l)7^0. 



(A.11) 



A. 2 The propagator H(x, /i) 
Let us consider the expression 



- ^(^°) 37^ ^(^') = - 1™ 4 ^(a:^ + A) 



(A.12) 



which is not a distribution in a four-dimensional space. To see this we consider the 
effect of integrating an arbitrary tempered test function with such an expression, 



/ d^x - Six"") ct>{x) ^ — [d^nr-4>l-—l d'n Hdr ^ 

J TT ^ ' d{x^) ^ '^^ ' At: J Jo r ^'*=" iir J Jo dt 



t=r 

(A. 13) 



where r = |a:|. The integral in the first term is, in general, divergent. The Hadamard 
finite part of the divergent integral drr~^ ip{r), where ip{r) is a tempered test 
function, is 



J^l — Loir) = lim / — Loir) + Ine o^fO) 
Jo r e-^0+ Je r 



°° dr 



(A. 14) 
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Taking (p{r) = J (PVt 01^^^, one can define the finite part of the integral (|A.13|) . In this 
way we can define a distribution, which we call T-'f 9{x^) ^7^21 



as 



J d^x Vf 



^Tpfd^x-eix''' 

J TT 



(t){x) 
d 



d{x'^ 



S{x'^) 0(x) 









L - 


r 


= lim 


f 


d"" 




J\x\>e 




= lim / 


'd^X 


"1 


e^0+ J 




TT 



1 

Air 



dt 



t=r 



1 
TT 



d{x'^ 



d 



\x\-e)-—^6ix^) + \ne5\x) 



0(x), (A.15) 



from which we see that Vf ^ 9{x^) 5{x'^) 



i^(x°)0(|x|-e)^5(x2) + 



Ine 6^(x) 



From eqs. ( |3.6|) and ( p.7|) we see that this distribution satisfies 



XaVf 



TT 



d 



6{x' 



TT 



9ix^) 



5{x'^) = Xa H(x, jj). 



(A.16) 



Notice that ( |A.16| ) is an equality between distributions: in the second member, the 
factor Xa kills the divergences. This equality shows that the propagator H{x, fl) must 
be equal to the distribution Vf ^ 6{x^) plus a distribution H{x, ft) which 

satisfies Xa H{x, Jj) = 0. It is easy to prove that such a distribution can only be a 
constant times the delta distribution, so we have 



H{x,fi)=Vf 



TT 



8[x^ 



(A.17) 



where C(/i) is a constant involving the arbitrary mass scale [i. Note that we could 
not have applied this argument to equation (|3^ ) because \ 9{x^) 5{x'^) is not 
a distribution. The constant C(/i) can be computed by evaluation of the integral 
/ d'^x H{x, Ji) 9{R — \x\), with R constant, using, on the one hand, expression ( |A.17| ) 
and, on the other hand, the following Fourier transform representation 



H{x, ji) 



d'^p 



-tpx 



In 



in 6{p^) sign{p^) 



(A.18) 



This can be easily shown to be equivalent to the definition of (x, /x) in eq. (|3.3| ) 
using lime^o+ In + ^c)] = In | | —iTx9{p^) and 2 6{—p^) — 1 = — sign(p°). We 
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obtain C(/x) =ln/i+7 — 1, so we finally have 
H{x,fi) = Vf 



+ [ln/x+7-1] 5\x), 
d 



lim <^ - ^(x") e{\x\ -e) — — 5{x') + [ln/ie+7-1] ^ . (A.19) 



This expression can be shown to be equivalent to 

■1 



H(x,fi)= \im\-9(x'')^6(x^ + \) + 

A^O- \ IT dX 



Inf-A/i" 1+7-1 



S\x) 



(A.20) 



which is, in fact, the definition for the propagator H used by Horowitz in ref. [22 



A. 3 Calculation of the integral I{x, y) 
We show here the calculation of the integral 

I(x, y) = d'^x' H{—x', Ji) 6{x + x')5{y + y') 



{A.21] 



which appeared in section 3. One can calculate this integral in two ways: using the 
representation ( |A.18|) for the propagator if as a Fourier transform or using expression 
( |A.lij| ). With the representation ( |A.l^ ) for H we have 



I{x,y) 



d'^p 



In 



P I 
— ■> 



(A.22) 



where p={p^,p^), x = {x, y) are here vectors in a two-dimensional space. That is, I is 
the Fourier transform of a logarithm in a two-dimensional space. The result, which 
can be found in ref. [0511, is 



Iix,y) = -Lrf(\) + (\n^ + ^]6ix)6iy) 



2lT 

— lim < 6(r — e) 



\n^ + ^)5ix)5{y) 



(A.23) 



where r = v^x^ + y'^. 

For completeness let us check that we can get this result using the representation 
( |Al9| ) for H, 



r 1 1 

/(x,i/) = £m |(ln/ie + 7-l)5(x)5(l/) + — dz ^^,^^,^3/, 



(A.24) 
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Introducing l = 6{r — e) + 6{e—r) in the second term one has 

f 1 1 

I{x, y) = hm j (ln/ie+7-1) 5{x)5{y) + — e{r-t) dz ^^^2)3/2 




hm |— - e{r-t) + (ln/ie+7-1) 5{x)5{y) + (l-ln2) ^^(x, y- e)| 



(A.25) 



where 



1 1 1 



27r l-ln2) r 



^(e-r) 1-^^ . (A.26) 



Now to prove (|A.23|) it remains to be seen that \mi^^Q+ Q{x,y]e) = 6{x)6{y). For 



this we use the following theorem (see, for example, ref. Suppose we have a 

function /(x; e), with x G iR", that satisfies the following conditions: 1) /(x; e) > 
for ||x|| < K, with K > some fixed constant; 2) in all the sets a < \\x\\ < ^, being 
a > an arbitrary finite constant, f{x;e) converges uniformly to zero as e ^ 0^; 
and 3) lim£^o+ /||x||<a, a>t /(^i ^) — 1- Then it can be shown that lime^o+ fi^] = 
S"'{x). On the one hand we see that Q{x,y; e)>0 and that Q{x,y] e) = for r>a>e, 
so, in the limit e — 0"*", Q{x, y; e) converges uniformly to zero for r > a > 0. On the 
other hand it is easy to prove that X-<^ ^^^dxdy Q{x,y;e) = 1, thus we have that 
lim,_o+ g{x,y;e)=5{x)S{y). 



Appendix B 

Riemann tensor in the linear approximation for /i^i, with cylindrical symmetry. 

It is easy to check that the Riemann components ( |2.14| ) for hf^^, with cylindrical 
symmetry can be written in terms of hi = hi, /i2 = hi, h^ = h^ + hg, h^ = h'[. — h^Q and 
h5 = hl = —hl in the following form: 

Rtztz = - h2,tt, R 'xzxz = - fcOS^6' /l2,rr + sin^6' j ) 



h 



2,r 



Rxzyz — ^ sin 2^ yl2,rr —J ; 

Rtxtx = 7 I ^3 « + COS 26* hi tt — "^ cos^ 6 — 2sin^6' — 4cos^6' /15 1^ — 4 sin^6' — ^ 
4 \ ' ' ' r ' r 
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Rtxty — ^ sin 29 yi4:,tt — ^l,rr H — ^ — 2 h^^tr H 
Rtxxy = T sin 6* i—hs^tr + ^4,tr + 2 ) , R 





R 



tzxz 



- COS 6 h 
2 



) 



Ryzyzi Rtytyi Rtzyz and -Rtj^yo; Can be obtained from the previous expressions inter- 
changing cos^ by sin^ in Rxzxz, Rtxtx, Rtzxz and Rtxxy respectively. 
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